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Coherent  eleetromagnelic  radiaiion  is  generated  by  the  oscillation  of  longitudinal  modes  associated  with 
optical  phonons  and  bulk  plasmons  in  InAs.  Dramatic  changes  to  the  radiation  spectra  occur  when  the  plasmon 
and  phonon  frequencies  are  in  close  spectral  proximity,  i.e.,  when  the  modes  become  strongly  coupled. 
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Plasmons  and  phonons  represent  fundamental  collective 
excitations  in  solids.  In  polar  materials,  the  plasmons.  and 
longitudinal  optical  (LO)  phonons  can  couple  to  form  hybrid 
modes  at  new  frequencies  that  were  first  observed  by  inco¬ 
herent  Raman  scattering  of  continuous  laser  light  from 
GaAs.1  The  availability  of  laser  pulses  with  duration  much 
shorter  than  the  oscillation  period  of  these  modes  allows  im¬ 
pulsive  excitation  of  coherent,  macroscopic  motion  of  charge 
carriers  and  lattice  ions  that  can  radiate  into  free  space." 
Emission  of  THz  frequency  electromagnetic  radiation  by  co¬ 
herent  phonons'-5*  and  bulk  plasmons'*  has  been  reported,  but 
only  limited  information  exists  for  die  situation  where  both 
modes  are  radiating.  Recent  experiments  have  provided  evi¬ 
dence  for  THz  emission  from  plasmon-phonon  modes  in 
InSb.  InAs.  and  GuAs.s  but  an  investigation  of  how  mode 
hybridization  affects  the  radiation  spectrum  has.  until  now. 
not  been  available. 

We  present  an  experimental  study  of  die  frec-space  radia¬ 
tion  originating  from  the  longitudinal  modes  of  a  polar  semi¬ 
conductor  as  a  function  of  coupling.  We  find  dial  the  THz 
frequency  emission  is  governed  by  bulk  plasma  oscillations. 
Only  when  the  longitudinal  phonon  vibrations  approach 
resonance  with  the  plasma  frequency — the  strong  coupling 
regime — do  they  significantly  affect  the  radiaiion  spectrum. 
In  addition,  we  observe  that  cold  plasmons  and  phonons  can 
rudiatc  coherently  even  w  hen  interacting  with  an  inhomoge¬ 
neous.  hot  electron-hole  plusmu  ut  a  density  several  orders  of 
magnitude  higher  than  the  background  doping  density.  Our 
experimental  results  arc  reproduced  by  model  calculations 
and  confirm  theoretical  predictions  of  Kuznetsov  and 
Stanton. ” 

Time -resolved  measurements  of  coupled,  coherent  oscil¬ 
lations  in  polar  semiconductors  have  been  primarily  made  by 
reflective  electro-optic  sampling  (REOS).1"  These  experi¬ 
ments  are  sensitive  to  electric  field  changes  that  occur  in  the 
narrow  (<1  /im)  surface  space-charge  region.  An  ul- 
trushori  laser  pulse  injects  electron-hole  pairs,  which  par¬ 
tially  screen  the  surface  field.  Long-range  Coulomb  coupling 
of  the  ultrafast  field  transient  to  both  the  electronic  polariza¬ 
tion  and  the  polar  lattice  can  initiate  collective  oscillations, 
which  weakly  modulate  the  surface  field  as  observed  in  re- 
llection.  A  complication  of  REOS  is  that  the  plasma  fre¬ 
quency  is  determined  by  the  sum  of  the  photocarrier  and 
background  carrier  densities.  Furthermore.  REOS  detects  os¬ 
cillations  that  are  confined  in  die  surface  field  layer;  Uiis  has 


been  shown  to  limit  the  minimum  plasmon  wave  vector  in 
InAs.11  Another  approach  is  to  measure  the  electromagnetic 
rudialion  emitted  by  coherent  macroscopic  dipoles  in  the 
bulk.  Two  popular  detection  techniques  are  optically  gated 
antennus  and  frec-space  electro-optic  sampling,  which  give 
the  temporal  phase  of  the  THz  pulse  relative  to  die  excitation 
pulse.  Phase  data  can  be  helpful  for  elucidating  the  starting 
mechanism  of  die  coherent  oscillations,  but  these  methods 
often  exhibit  dispersion  and  Reststrahl  absorption  at  optical 
phonon  frequencies  in  1I1-V  semiconductors. 121 '  Here,  we 
use  a  phase-insensitive  arrangement  to  make  broadband  mea¬ 
surements  of  die  emitted  radiation  spectrum,  which  elimi¬ 
nates  signal  distortion  problems  in  the  spectral  region  of  in¬ 
terest  between  1-10  THz. 

Our  experiments  are  based  on  an  interferometric  tech¬ 
nique  used  by  Kersting  et  al.  to  study  coherent  plasmon  os¬ 
cillations  in  bulk  GaAs.1’  Two  30  fs  pump  pulses  derived 
from  a  Tirsapphire  laser  (center  wavelength  770  nm)  are  fo¬ 
cused  on  slightly  separated  spoLs  on  an  InAs  crystal  (£k> 
=  0.36  eV  at  300  K)  at  an  incidence  angle  of  45°.  Spatial 
separation  of  the  pump  pulses  represents  an  important  modi¬ 
fication  of  the  original  setup;  each  excited  spot  generates  an 
independent  THz  frequency  pulse  without  complications 
from  pump  beam  overlap.11  The  THz  pulses  propagate  into 
free  space  and  are  collected  by  an  off-axis  parabolic  mirror, 
collimated,  and  spatially  filtered.  A  second  mirror  focuses  the 
two  THz  beams  on  a  liquid-He  cooled  Si  bolometer.  By  ad¬ 
justing  the  relative  time-delay  of  die  two  pump  pulses,  the 
detector  signul  maps  out  a  linear  inlcrfcrogram  of  the  elec¬ 
tromagnetic  transient  emitted  from  the  semiconductor.  Al¬ 
though  the  THz  pulses  originate  from  spatially  separate  re¬ 
gions.  their  diffraction  angle  is  sufficiently  lurge  to  allow 
interference  in  the  fur  field.  The  pump  irradiancc  is  kept  low' 
to  suppress  generation  of  fur-infrared  light  by  surface 
rectification.15  The  InAs  samples  are  of  two  types:  (i)  0.5- 
mm  and  1 -nun-thick  bulk  LEG  grown  wafers  oriented  (111) 
or  (1(H))  and  doped  with  S  or  Tc  and  (ii)  ~  1  -/zm- thick  MBE 
layers  doped  with  Si  and  grown  on  intrinsic  (100)  InAs  sub¬ 
strates.  Other  than  doping  concentration,  we  could  not  iden¬ 
tify  any  dependence  on  sample  orientation  or  preparation 
method  in  the  THz  data. 

Figure  1  displays  (a)  time-resolved  linear  interferograms 
and  (b)  corresponding  Fourier- transform  spectra  obtained  ut 
room  temperature  with  samples  al  lour  progressively  higher 
donor  densities  while  maintaining  the  same  excitation  irradi- 
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FIG.  1.  (a)  Time-resolved  linear  interferograms  for  bulk  (1(H)) 
//-doped  InAs  samples  and  (h)  corresponding  Fourier  transform 
spectra  (solid  lines).  Labeled  doping  densities  are  in  units  of  cm  \ 
The  dashed  vertical  line  indicates  the  bare  LO  phonon  frequency. 
Doited  curves  in  (b)  are  calculated  using  Eq.  (2). 

ancc.  The  lop  spectrum  with  a  lightly  //-doped  hulk  wafer 
reveals  a  broad  plasma  peak  centered  at  i'/,=2.1  THz  and  u 
much  weaker  peak  near  the  bare  LO  phonon  frequency  of 
r'l  o=7.3  THz  (dashed  vertical  line  lor  reference).  As  the 
donor  density  increases,  three  features  are  evident:  (i)  the 
plusma  peak  moves  to  higher  frequency,  (ii)  coupling  of  the 
plusmon  like  to  the  phononlike  mode  strengthens  causing  a 
blucshifl  of  the  upper-hybrid  mode,  und  (iii)  the  spectral 
weight  of  the  upper-hybrid  mode  dramatically  increases  ac¬ 
companied  by  signilicanl  broadening.  The  bottom  spectrum, 
where  the  sample  is  doped  to  be  very  close  to  the  antieross- 
ing  frequency  (i.c..  where  v,,=  vLO).  exhibits  a  pronounced 
shift  of  the  upper-hybrid  mode.  The  plasma  frequency  is  just 
above  the  bare  LO  phonon  frequency  in  this  sample  so  we 
identify  the  upper-hybrid  mode  as  the  plasmonlike  mode. 

The  spectrul  locations  of  both  the  upper-  and  lower-hybrid 
modes  are  explained  by  the  concentration  of  ionized  donor 
atoms  in  the  bulk  sample.  The  longitudinal  mode  frequencies 
of  the  coupled  plusmon-phonon  system  are  zeros  of  the  di¬ 
electric  response  function  in  the  weak  damping,  long  wave¬ 
length  limit,  i.c..  where 


i-^i  + 


c0— e, 


0)~  \—w~/wfv~ 


=  0. 


(I) 


The  low-  and  high-frequency  dielectric  constants  arc  e„ 
=  11.8  and  ew=  14.55.  respectively.  wlu/27r=6.54  THz  is 
the  transverse  optical  phonon  frequency,  and  the  plasma  fre¬ 
quency  is  <i>~= 4  7TiVe'/ //?*£„ .  The  optical  effective  mass 
(wr*)  accounts  for  the  nonparabolic  conduction  band  of 
InAs.  which  becomes  more  important  as  the  doping  density 
(AO  increases.  1 1  The  mass  of  the  donor  ions  is  taken  lo  be 
inlinile.  Results  of  this  calculation  are  shown  in  Fig.  2  with 
solid  lines.  The  filled  circles  are  obtained  from  the  THz  spec¬ 
tra  of  six  different  //-doped  samples. 


FIG.  2.  Dispersion  of  the  plasmon-phonon  modes  of  InAs  as  a 
function  of  donor  doping.  Filled  circles  are  front  experiment;  solid 
curves  are  derived  from  Eq.  (1). 

The  center  frequencies  of  the  spectral  peaks  do  not 
change  as  the  optical  excitation  power  is  varied  bv  a  factor  of 
3.  showing  the  coherent  oscillations  are  independent  of  the 
level  of  photocarrier  injection.  This  is  in  agreement  with  pre¬ 
viously  reported  results  on  THz  emission  by  cold  plasmons 
in  GaAs  (Ref.  6)  und  InSb.  The  excitation  pulses  are  p  po¬ 
larized  and  the  detected  signuls  exhibit  no  dependence  on 
uzimulhul  rotation  for  either  (100)  or  (111)  InAs  samples. 
This  indicates  the  oscillutions  are  initialed  by  ultrufusl 
screening  of  the  surface  field  by  optically  injected  carriers,  u 
mechanism  that  is  well  known  for  GaAs."1  Recent  REOS 
measurements  of  coherent  phonons  in  InAs  are  also  consis¬ 
tent  with  displacive  excitation  by  surface  field  screening.11 

The  absorption  depth  of  the  pump  light  is  >200  nm. 
which  causes  uniform  screening  of  the  surface  accumulation 
field  (depth  **5  nm.  Ref.  17)  und  also  penetrates  into  the 
neutral  bulk  region,  i.c..  the  region  where  the  coherent  oscil¬ 
lations  originate.  We  estimate  the  optically  excited  carrier 
density  in  the  pump  absorption  region  to  be  well  in  excess  of 
I0IX  cm-  '.  This  is  at  least  an  order  of  magnitude  greater 
density  than  the  background  doping  concentration,  yet  the 
optically  injected  electron-hole  pairs  do  not  affect  the  fre¬ 
quency  of  die  cold  plasma  oscillations  in  our  THz  spectra. 
Absence  of  a  coherent  contribution  from  the  photocarriers  is 
due  to  their  inhomogeneous  (i.c..  quasiexponential)  distribu¬ 
tion  in  the  bulk  region  adjacent  to  die  surface. IS 

A  quantitative  understanding  of  the  THz  signals  is  ob¬ 
tained  with  the  time-domain  model  of  Ref.  9.  which  couples 
the  electronic  polarization  with  die  lattice  displacement.  Fou¬ 
rier  analysis  gives  the  power  spectrum  of  die  emitted  radia¬ 
tion 

F(w)o:to,|£'(<«»)|2 

a»;(e*D  |  ( tu)  —  4  7rH )  +  fi[  D  .(w)  —  w"J 
e  XD ,  (  io  )  D  2(  to )  —  4  7rft  a>“ 
Dl(o>)=-io1  +  iyrhio+o)lo:, 

Di(o>)  =  —  fti2+/ to~  (2) 
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with  wLO=  27ri'LO  and  ft  =  wTO-(e„— ew)/4ir.  The  electric 
held  transient  at  the  edge  of  the  space-charge  layer  is  mod¬ 
eled  as  £(/)  =  £„[tanh(//r)+ 1],  where  r  uccounls  for  the  re¬ 
sponse  time  of  the  surface  field  screening.  This  function  is  an 
analytic  approximation  to  the  calculated  surface  field 
dynamics.61119  The  Fourier  transform  of  £(/)  gives  the 
spectrum  of  the  surface  transient  |£(oj)|2«csch2(ira>T/2), 
for  w>0.  The  width  of  the  accumulation  layer  scales  ap¬ 
proximately  as  \/\[N  so  we  expect  luster  screening  as  the 
donor  density  increases.20  The  bare  phonon  relaxation  time  is 
set  at  \/y !,/,  =  !  ps  reflecting  its  anharmonic  decay1  and  best 
lils  to  the  experimental  data  are  obtained  with  a  density- 
independent  plusmon  damping  lime  (i.e..  electron  momen¬ 
tum  relaxation  lime)  of  l/yt|=  125  Is.  Results  of  the  model 
calculations  are  depicted  by  dotted  curves  in  Fig.  1(b).  The 
surface  field  screening  response  times  are  found  by  lilting  the 
amplitude  of  the  upper-hybrid  mode  peaks,  yielding  values 
of  t—  76.  68.  58.  and  52  Is  for  the  increasing  sequence  of 
doping  densities. 

In  the  strong  coupling  regime,  the  phononlike  mode  lakes 
on  the  damping  rate  of  the  plasmons  while  the  plasmons 
become  more  phononlike.22  Most  experiments — in  particular 
time-resolved  Raman  measurements — study  the  behavior  of 
just  one  mode  and  deduce  mode  hybridization  indirectly.2' 
Here,  we  can  clearly  resolve  the  spectra  of  both  modes  as  a 
function  of  coupling  and  directly  address  this  question.  It  is 
seen  that  the  contribution  of  the  phononlike  mode  to  the 
radiated  THz  signal  depends  critically  on  the  spectral  prox¬ 
imity  of  the  plasma  frequency.  When  the  modes  are  strongly 
coupled,  the  lattice  dynamics  can  significantly  affect  the 
electric  Held  that  is  responsible  for  the  emitted  radiation.  We 
believe  our  data  represents  the  first  observation  of  this  effect, 
which  was  predicted  in  Ref.  9. 

We  now  discuss  the  spectral  linewidlhs  in  Fig.  1(b).  The 
decay  of  weakly  coupled  coherent  plasmons  is  associated 
with  the  electron  momentum  relaxation  time  and  therefore 
the  mobility  of  the  semiconductor.  As  with  GaAs.  the  mobil¬ 
ity  of  room  temperature  InAs  should  be  dominated  by  LO 
phonon  scattering.1 1  The  plasmon  dumping  rate  derived  from 
our  experiments  is,  however,  uboul  3  times  higher  than  indi¬ 
cated  by  Hall  mobility  measurements.  THz  inter ferograms  at 
80  K  (not  shown)  reveal  that  this  unexpectedly  high  plasmon 
damping  rale  does  not  change  when  compared  to  room  tem¬ 
perature  spectra.  Because  LO  phonon  scattering  is  strongly 
temperature  dependent,  it  must  therefore  play  u  minor  role 
here.  A  constant  plasmon  damping  rate  is  used  in  the  model 
and  lils  the  spectra  well  at  all  doping  densities,  which  shows 
that  ionized  impurity  scattering  is  negligible.  Nonuniform 
spatial  excitation  of  the  coherent  volume  due  to  the  trans¬ 
verse  laser  profile  can  cause  inhomogeneous  broadening.9'24 
but  we  found  the  plasmon  line  shape  to  be  insensitive  to 
spatial  filtering  of  the  THz  beams.  Landau  damping  can  be 
ruled  out  because  large  wave  vector  modes  do  not  efficiently 
radiate.14  If  the  coherent  plasmons  are  confined  in  the  surface 
accumulation  region  (</^2Xl()'’  cm-1)  as  observed  in 
RHOS.  single  particle  excitation  could  lake  pluce  and  the 
plasmons  would  be  severely  broadened  by  Landau  damping 
accompanied  by  strong  emission  at  the  bare  LO  phonon 
frequency.16-2'  We  find  no  peaks  at  the  LO  phonon  frequency 


FIG.  3.  Line  shape  of  the  plasma  peak  at  two  different  excitation 
levels.  The  solid  curve  is  obtained  at  approximately  twice  the  irra- 
diance  of  the  dashed  curve. 

in  the  THz  spectra  of  the  samples  n-doped  at  1 .8  X  101  and 
3.1  X101  cm”'.  It  can  be  concluded  that  longitudinal 
modes  with  wave  vector  </**()  are  emitted,  with  a  corre¬ 
sponding  plasmon  displacement  amplitude  of  at  least  50  nm. 
The  long- wavelength  nature  of  the  radiating  modes  means 
that  Landau  damping  of  the  coherent  plasmons  is  not  rel¬ 
evant. 

The  high  plasmon  damping  rale  is  explained  by  the  pres¬ 
ence  of  hot  carriers  1.2  eV)  injected  by  the 

near-infrared  pump  pulse.  Although  the  inhomogeneous  lon¬ 
gitudinal  distribution  of  the  photocarriers  does  not  contribute 
to  the  coherent  signals,  previous  work  has  revealed  that  there 
can  be  significant  damping  of  coherent  plasmons  and 
phonons  due  to  scattering  with  optically  generated  electrons 
and  holes.  l8-22~24-26  To  confirm  the  presence  of  strong  scat¬ 
tering  by  nonequilibrium  carriers,  the  excitation  irradiancc  of 
the  near-infrared  pump  light  is  increased  by  about  a  factor  of 
2.  leading  to  a  larger  photocarrier  density  in  the  radiating 
volume  beneath  the  surface  accumulation  layer.  Changes  in 
the  THz  spectrum  in  the  vicinity  of  the  coherent  plusmon  for 
the  sample  doped  6.3X  l()lb  cm-3  are  depicted  in  Fig.  3. 
With  higher  photocarricr  density  (solid  line),  we  observe  in¬ 
homogeneous  broadening  on  the  high-frequency  side  of  the 
plasmon  peak.  Moreover,  the  model  described  by  Eq.  (2) 
assumes  homogeneous  broadening  and  thus  the  calculated 
curves  consistently  underestimate  the  experimental  spectra 
on  the  high-frequency  wings  of  the  peaks  in  Fig.  1(b).  This 
strongly  indicates  inhomogeneous  broadening  by  optically 
injected  electron-hole  pairs  is  an  important  factor  in  the  in¬ 
terpretation  of  our  data.  The  extracted  plasmon  damping  time 
of  l/yci~125  fs  is  also  consistent  with  cold  electrons  scat¬ 
tering  with  holes  of  a  much  higher  density.27  Note  thul  the 
hot  electron-hole  distribution  induced  by  abovc-bandgap  la¬ 
ser  pulses  persists  for  3-4  ps  in  InAs.  i.e..  longer  than  the 
time  scale  of  our  measurements.2' 

In  summary,  we  have  studied  the  emission  of  far-infrared 
electromagnetic  transients  from  bulk  InAs  with  a  high  fre¬ 
quency.  uniform  bandwidth  detection  scheme.  Distinguish- 
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able  contributions  from  both  coherent  plasmons  and  phonons 
are  observed  in  the  radiated  signals.  By  varying  the  donor 
doping  and  hence  the  plasma  frequency,  the  long  range  Cou¬ 
lomb  coupling  of  the  electronic  polarization  with  the  polar 
lattice  changes,  resulting  in  striking  variation  of  the  emission 
spectra.  A  nonequilibrium  electron-hole  plasma  injected  by 


the  excitation  laser  pulse  is  primarily  responsible  for  damp¬ 
ing  the  coherent  oscillations. 
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